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Abstract. It is known that Struve function H„ and modified Struve function L„ are closely connected 
to Bessel function of the first kind J„ and to modified Bessel function of the first kind /„ and possess 
representations through higher transcendental functions like generalized hypergeometric 1F2 and Meijer 
G function. Also, the NIST project and Wolfram formula collection contain a set of Kapteyn type 
series expansions for li v (x). In this paper firstly, we obtain various another type integral representation 
formulae for L„(:r) using the technique developed by D. Jankov and the authors. Secondly, we present 
some summation results for different kind of Neumann, Kapteyn and Schlomilch series built by I v (x) 
and L„(i) which are connected by a Sonin-Gubler formula, and by the associated modified Struve 
differential equation. Finally, solving a Fredholm type convolutional integral equation of the first kind, 
Bromwich- Wagner line integral expressions are derived for the Bessel function of the first kind J„ and 
for an associated generalized Schlomilch series. 

1. Introduction 

Bessel and modified Bessel function of the first kind, Struve and modified Struve function possess 
power series representation of the form [63] : 



U*) - E krzfhurv W = E 



n>Q y ' n>0 y ' 

I2J 



= E r ^ ( l„ X 3 V m*) - E 



^ o r(n + |)r(n + 2 , + |)' Z-.r( n + l)r(n + v 



where v e M. and z G C. Struve [56] introduced H„ function as the series solution of the nonhomoge- 
neous second order Bessel type differential equation (which carries his name). However, the modified 
Struve function L„ appeared into mathematical literature by Nicholson [40, p. 218]. Applications of 
Struve functions are manyfold and include among others optical investigations [60, pp. 392-395]; gen- 
eral expression of the power carried by a transverse magnetic or electric beam, is given in terms of 
L n _|_i [2]; triplet phase shifts of the scattering by the singular nucleon-nucleon potentials oc exp(— x)/x n 
[20]; leakage inductance in transformer windings [26]; boundary element solutions of the two-dimensional 
multi-energy-group neutron diffusion equation which governs the neutronic phenomena in nuclear reactors 
[27]; effective isotropic potential for a pair of dipoles [36]; perturbation approximations of lee-waves in a 
stratified flow [37]; quantum-statistical distribution functions of a hard-sphere system [42]; scattering of 
plane waves by circular cylinders for the general case of oblique incidence and for both real and complex 
values of particle refractive index [54]; aerodynamic sensitivities for subsonic, sonic, and supersonic un- 
steady, non- planar lifting-surface theory [14]; stress concentration around broken filaments [19] and lift 
and downwash distributions of oscillating wings in subsonic and supersonic flow [61, 62]. 

Series of Bessel and/or Struve functions in which summation indices appear in the order of the consid- 
ered function and/or twist arguments of the constituting functions, can be unified in a double lacunary 
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form: 

n>l 

where x ^ £j(x) = \ij + ajx, j e {1, 2}, x e {0, 1, . . . }, z e C and 8S U is one of the functions J v , I Vl H u 
and L„. The classical theory of the Fourier-Bessel series of the first type is based on the case when 
8§ u = J„, see the celebrated monograph by Watson [63]. However, varying the coefficients of £\ and £2, 
we get three different cases which have not only deep roles in describing physical models and have physical 
interpretations in numerous topics of natural sciences and technology but are also of deep mathematical 
interest, like e.g. zero function series [63]. Hence we differ: Neumann series (when a\ 7^ 0, a 2 = 0), 
Kapteyn series (when a\ ■ a 2 7^ 0) and Schlomilch series (when a\ = 0,a 2 7^ 0). Here, all three series 
are of the first type (the series' terms contain only one constituting function 8$„); the second type series 
contain product terms of two (or more) members - not necessarily different ones - from J„, and 
~L U . We also point out that the Neumann series (of the first type) of Bessel function of the second kind 
Y v , modified Bessel function of the second kind K v and Hankel functions (Bessel functions of the third 
kind) H^\H^ have been studied by Baricz, Jankov and Pogany [5], while Neumann series of the second 
type were considered by Baricz and Pogany in somewhat different purposes in [6, 7]; see also [31]. An 
important role has throughout of this paper the Sonin-Gubler formula which connects modified Bessel 
function of the first kind I v , modified Struve function L v and a definite integral of the Bessel function of 
the first kind J„ [23, p. 424] (actually a special case of a Sonin-formula [63, p. 434]): 

it -n f°° J "( ax ) dx n IJ I \ T I w 

where Sft(v) > — 2 , a > and Sft(n) > 0; see [63, p. 426] for the historical background of (1.1). 
Thus, under extended Neumann series (of Bessel J„ see [63]) we mean the following 

^O*) ~J20n8§„ n+ri (ax), 

n>l 

where 8S V is one of the functions I v and L„ . Integral representation discussions began very recently with 
the introductory article by Pogany and Siili [47], which gives an exhaustive references list concerning 
physical applications too; see also [5]. In Section 2 we will concentrate to the Neumann series 

(1-2) ^n,v( x ) : = Pnl^n+ n {ax) . 

Secondly, Kapteyn series of the first type [32, 33, 41] are of the form 

n>l 

more details about Kapteyn and Kapteyn-type series for Bessel function can be found also in [4, 6, 15, 57] 
and the references therein. Here we will consider specific Kapteyn-type series of the following form: 

(1.3) ■■= E ^ _ L -( xn )) ; 

n>l 

this series appear as auxiliary expression in the fourth section of the article. Thanks to Sonin-Gubler 
formula (1.1) we give an alternative proof for integral representation of &"^(x), see Section 3. Thirdly, 
under Schlomilch series [52, pp. 155-158] (Schlomilch considered only cases e {0,1}), we understand 
the functions series 

(1.4) 6^(z):=5>n^((^ + ^)- 

n>l 

Integral representation are recently obtained for this series in [30], summations are given in [59]. Our 
attention is focused currently to 

6 £'» ^J2^(M x n)-U(xn)) . 

n>l 
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The next generalization is suggested by the theory of Fourier series, and the functions which naturally 
come under consideration instead of the classical sine and cosine, are the Bessel functions of the first kind 
and Struve's functions. The next type series considered here we call generalized Schlomilch series [63, p. 
622], [28, p. 1803] 

&o ( x \~ v a n J v (nx) + 6„H ly (nx) 



2T(v + 1) 

V ; n>l 

For further subsequent generalizations consult e.g. Bondarenko's recent article [9] and the references 
therein and Miller's multidimensional expansion [38]. A set of summation formulae of Schlomilch series 
for Bessel function of the first kind can be found in the literature, such as the Nielsen formula [63, p. 
636]; further, we have [58, p. 65], also consult [21, 45, 50, 59, 64, 65]. Similar summations, for Schlomilch 
series of Struve function, have been given by Miller [39], consult [59] too. 

Further, we are interested in a specific variant of generalized Schlomilch series in which J„,H„ are 
exchanged by I v and L„ respectively, when a n ,b n arc of the form a = 0, a n = 2^ v n v ^ ti x v = —b n , 
H > v > 0, which results in 

Zi,L( x ) .= V" Iyjnx) - L v (nx) 

n>l 

Its alternating variant T^(x) we perform setting (— l)™ _1 a„ i->- a n , where n£ {0, 1, . . . }: 

■= E {IAnx) u(nx)) ■ 

Summations of these series are one of tools in obtaining explicit expressions for integrals containing 
Butzer-Flocke-Hauss complete Omega-function fi(x) [10, 11, 12] and Mathieu series S(x),S(x) [35, 49]. 

Rayleigh [51] has showed that series &q v {z) play important roles in physics, because they are useful 
in investigation of a periodic transverse vibrations uniformly distributed in direction through the two 
dimensions of the membrane. Also, Schlomilch series present various features of purely mathematical 
interest and it is remarkable that a null-function can be represented by such series in which the coefficients 
are not all zero [63, p. 622]. Summation results in form of a double definite integral representation for 
&^ v {z), achieved via Kapteyn-series, have been recently derived in [29]. 

Finally, we mention that except the Sonin-Gubler formula (1.1) another main tool we refer to is the 
Cahen formula on the Laplace integral representation of Dirichlet series. Namely, the Dirichlet series 



n>l 



where 3?(r) > 0, having positive monotone increasing divergent to infinity sequence (A n ) n >i, possesses 
Cahen's integral representation formula [13, p. 97] 

f°o roo /■[A _1 (t)] 

(1.6) $ a ( r )=r e~ rt V a n dt = r / <3 u a(u)dtdu, 

where d x := 1 + {x}^. Here, [x] and {x} = x — [x] denote the integer and fractional part of x e M, 
respectively. Indeed, the so-called counting sum 

n: \„<t 

we find by the Euler-Maclaurin summation formula, following the procedure developed by the second 
author [46], see also [49]. Namely 

[A _1 (*)] A\-Ht)] 
*f a (t) = V a n = / Q u a(u)du, 
n=i Jo 

since A: R + M + is monotone, there exists unique inverse A -1 for the function A: M + M + , A|n = (A n ). 
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2. L„ as a Neumann series of modified Bessel I functions 

Let us observe the well known formulae [44, Eqs. 11.4.18-19-20] 

(2n + v+l)T(n + v + l) 



H„(*) 



v^r(,+ i) ^ n!(2n+l)(2n + 2,+ l) 



O-rr X! „W^\ 1\ J n+v+±{ z ) 



(§)" 



1\ "Ai+A 



where the first formula is valid for — v £ N. So, having in mind that Iv(z) = — i 1_ "H„(i,2) and J v {iz) — 
i v I v (z), we immediately conclude that 

4 vM-l) n (2n + z/+l)r(n + i/ + l) 



(2.1) 



L„(z) = < 



^O'+dS n!(2n + l)(2n + 2^+l) 



'2n+ 



*+l(*0 



2tt E 



n>0 



d („+§) 



1\ 



(I) 



(-§)" 



.r^ + i) ^o n! ( n + z/ + 5) 

However, all three series expansions we recognize as Neumann-series built by modified Bessel functions 
of the first kind. This kind of series have been intensively studied very recently by the authors and D. 
Jankov in [5]. Exploiting the appropriate findings of that article, we give new integral expressions for the 
modified Struve function L„. 

First let us modestly generalize [5, Theorem 2.1] which concerns OTi^x), to integral expression for 
O^.r, defined by (1.2), following the same procedure as in [5]. 

Theorem 1. Let (3 e C 1 (M + ), /3|n = {/3„} nG N; M > an d assume that 
(2.2) W ^ " 



lim 



< 



n— >oo 71 



Then, for fi, n such that minjr; + §, [i + r\ + 1} >0 and 



x€ ^0,2min ^1, {{e/^f limsupn^|/3 n | 1/n ^ := J & . 



we have the integral representation 



r (pa + v + 2) 



duds . 



(2.3) m^ix) = - ^ A-fr^u + n+^j d s ^ 

Proof. The proof is a copy of the proving procedure delivered for [5, Theorem 2.1]. The only exception 
is to refine the convergence condition upon 9l M)?J (a:). By the bound [3, p. 583]: 

:§r 



I> + 1) 



e 4(^+1). 



where x > and v + 1 > 0, we have 



i*w*)i < (f )" + " e 3 ^^ e 

n> 



jrfrm + fj + l) ' 
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so, the absolute convergence of the right hand side series suffices for the finiteness of Vi^^ix) on J* p. 
However, condition (2.2) ensures the absolute convergence by the Cauchy convergence criterion. 

The remaining part of the proof mimicks the one performed for [5, Theorem 2.1], having in mind that 
jj, = 1 reduces Theorem 2 to the ancestor result [5, Theorem 2.1]. □ 



Theorem 2. If v > and x € (0,2), then we have the integral representation 
2T(v + 2) 

V*Fr(i/ + i) 



Iv (x) ; 3 , I v +l{x) 



poo r[u] g / of \ \ 

(2.4) =^ ^ - ( r(2. + , + l) Wl ( s) ),^_ai_^jd«d fl) 

e i7rs (2s + i/+ l)r(.s + i/+ 1) 



where 



V^r (i/ + i) r( s + 1) ( s + 1) ( a + 1/ + 1) ■ 

Proof. Consider the first Neumann sum expansion of L„(x) in (2.1), that is 

t / \ _ 4 (-l)"(2n + I /+l)r(n + !/+!) 

1 J ~>/5Fr(i/+|);^ n!(2n+l)(2n + 2^ + l) ^+"+^ x > 

2i> + 2) 

■U+i(x) 



v^i> + §) 



4 x ^ {-l) n -\2n + v+ l)r(n + u + 1) r 



V^r(i/+^)^ n!(2n + l)(2n + 2z/+l) 



Observe that 



in which we specify 



Since 



2T{v + 2) 



(-l)"~ 1 (2n + v + l)r(n + ^ + 1) 
(j/ + |) r(n + 1) (n + |) (n + i/ + §) ' 



2.s 



i/-2 



we deduce (by means of Theorem 1) that (2.4) is valid for x £ J?p = (0, 2). □ 
Theorem 3. For v + 2 > and x <G (0, 2) we /icwe i/ie integral representation 

(2.5) - y^/^jW 

(2.6) =/°°/ H i(r(. + » + .)i, t< (.))..( f5 «_)«,, 

w/iere 

e-(f) s 



/3( S ) 



2^ r( s + i)( s + i) 

Proof. Let us observe now the second Neumann sum expansion of L„(x) in (2.1): 



n>l V 



x). 
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In other words, 



in which we specify 



L„(x) 



(f) f 



2ir r(s + l)(s+i) 

The convergence condition (2.2) reduces to the behavior of the auxiliary series 



gr(n + , + i) 



— 1F2 

TT 



1 

2' ^ ^ 2 



which is convergent for all bounded x € C, unconditionally upon Here 1F2 denotes the hypergeometric 
function defined by series [1, p. 62] 



1-F2 



a 

61, b 2 



= E 

n>0 



(a)n 



(6l)„(6 2 )„ 



However, for v > —2 we have the integral expression [63, p. 79] 



(2.7) 



V^r(z. + i) J 



(l-t 2 y-i cosh(zt)dt, 



where z e C and > — |. This was used in the proof of the ancestor result (2.3), see [5, Theorem 

2.1]. Now, we apply Theorem 1 and conclude that (2.5) is valid for (0, 2). □ 

The third formula in (2.1) one reduces to the case <i yi 1 i(x). However, we shall omit the proof, since 
the slightly repeating derivation procedure used for (2.5) directly gets the desired integral expression. 

Theorem 4. If v + 2 > and x e (0, 2), then we have the integral representation 



~L v {x) - 



x v sinhx 

2«Vir (!/ + §) 



r°° /■["] d / \ 
I I ^(r(« + DM(*))j. 



T(s + 1) 



du ds . 



where 



r(i/ + |) r(« + i) (« + !/ + 1) ' 

Now, applying the integral representation (2.7) we derive another integral expression for L v (x) in 
terms of hypergeometric functions in the integrand. 

Theorem 5. Let v > — h. Then for x > we have 



^F2^-lr(^ + |)r(| + l)r(f + f) J a 

x 4F5 



2' 2 > u ^ 2 > ^ 



3 ^+1 1/ _i_ 3 v_ 1 5 . . _i_ 3 
2' 2 '2'4'2'4''4 



(1 — 1 2 )" + 3 cosh(rf) 



where 



4-T5 



ai, a2, 03, a.4 
61,62,63,64,65 



E 



4 

II ( a j)« . 
3=L ^_ 

r " n! 



II (6j)n 

stands for the generalized hypergeometric function with four upper and five lower parameters. 
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Proof. Consider the first Bessel function series expansion for li„(x) given in (2.1). Applying mutatis 
mutandis the integral representation formula (2.7), the Pochhammer symbol technique, the familiar 
formula (A) n (n + A) = A(A + 1)„, n e {0, 1, . . . }, and the duplication formula 

r(2zH^Jr(z)r(z + i) 

to the summands, we get the chain of equivalent legitimate transformations: 

(-l)"(2n + i/ + l)r(n + i/+l) 2(f) 2 " 



^ 2n +"+3 cosh(art) d* 



(I-/- 



4(f) 

(n + ^) r(n + v + 1) [-£ (1 - t 2 



?> / (l - 1 2 )"+5 cosh(^) 
("+2) ^0 

E 



n>0 



(n + |) (n + z/ + |) T (2n + z/ + §) n\ 



■At 



4(f) I/+1 ( v + i)r(i/ + i) 



V5F(i/ + |)r(i/ + ^)r(f + f)r(| 

(5/)„(^)n(^+^ n ^+l)n 



E 



16 



/ (1 - cosh(:ri) 
Jo 

(l_^2l 



di, 



, „ (i) B m B (5+iU5+s) B ("+i) B »i 

which proves the assertion. 

By virtue of similar manipulations presented above, we conclude the following results. 
Theorem 6. Let v > — \ and x > 0. Then there holds 



□ 



L„(a;) 



-1 



2"- 1 7rI> + l) 7 

,L/+1 



.7,' 



(1 - t 2 ) v cosh(a;i) iF 2 
•1 / 



V5F2"+*r(i/ + §) 7 



cosh(xi) iF 2 



2 ' 
1 
2 



■(1-t 2 ) dt 



-(l-t 2 ) dt 



The proof of Theorem 6 follows from the same proving procedure as the previous theorem but now 
considering the second and third series expansion results in (2.1), so we shall omit the proofs of these 
integral representations. 

3. Integrals containing Q(x) —function and Mathieu series via summation of T 7 j, L (x) 

By virtue of the Sonin-Gubler formula (1.1) we establish the convergence conditions for the generalized 
Schlomilch series T^(x) and T^(x). As for n enough large we have 

2n v - 1 f°° J u {ax)dx 



(3.1) 



I v {an) — L„(an) 



(1 + n 2 x 2 )x' / 



'( n ) > 



we immediately conclude that the following equiconvergences hold true 

~ C(m - v + 1), %l%(x) ~ »?(a* - v + 1) , 

that is, ll'^(x) converges for fi > v > 0, while Tf/j^a;) converges for \i + 1 > t/ > 0. On the other hand, 
we connect Tf/J^x) and the Butzer-Flocke-Hauss (BFL) complete Omega function [10, Definition 7.1] 



Q(w) — 2 / sinh(ww) cot (701) du, 



w € 
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By the Hilbert transform terminology, fl(w) is the Hilbert transform J$?(e wx )i(0) at of the 1-periodic 
function (e~ wx )i defined by the periodic continuation of the following exponential function [10, p. 67]: 
e~ xw , x e [-§, §) , w e C, that is, 



^(e- xu, )i(0) :=P.V.y* 2 e"™ cot(7ru) du = 



where the integral is to be understood in the sense of Cauchy's Principal Value at zero, see e.g. [12, 48]. 

On the other side by differentiating once (1.1) with respect to n we get a tool to obtain Mathieu series 
S(x) (introduced by Emile Leonard Mathieu [35]) and its alternating variant S(x) (introduced by Pogany, 
Srivastava and Tomovski [49]), which are defined as follows 

s(x) y 2n s(x) - V 2 (~ 1 )"" lw 

{> ^(x 2 +n 2 ) 2 ' (> 2^( x 2 + n 2y 

n>l n>l 

Closed integral expression for S(r) was considered by Emersleben [18] and subsequently by Elbert [17], 
while for S^x) integral representation has been given by Pogany, Srivastava and Tomovski [49]: 

(3.3) *)=if^. 

x J e* + 1 

Another kind integral expressions for underlying Mathieu series can be found in [49]. 
Theorem 7. Assume that $t(i/) > and a > 0. Then we have 

r°° Max) n(27rx) dx ^ v r°° e _„ t [ [e \ u t el .u ( ^ Aau) _ Iv{au)) - ) dtdu _ 

Jo x v smh(7rx) Jo Jo 

Proof. When we multiply (1.1) by (—l) n ~ 1 n and sum up all three series with respect to n € N, the 
following partial-fraction representation of the Omega function [10, Theorem 1.3] 



■di, 



immediately give 

J u (ax) Q(2irx) 



I 



o x v sinh(7rx, 

n>l 



nn(2nw) _ ^ 2(-l)"- 1 n 
sinh(7rw) ^ n 2 + w 2 

x ' n=l 



Ax = Y J {-V) n ~ 1 n- v (I v (an) - L v (an)) = X#(a) . 



We recognize the right-hand-side sums as Dirichlet series of I v and L v , respectively. Being 
^2{-l) n - 1 n-'I v (an) = ^ e i7r ("- 1 )j i ,(an) e - I/ln ", > 0, 

n>l n>\ 

we get 

^(-l) n - 1 n- v I v {an) =v / e~ vt ^ e i7r( "- 1) / i ,(an) dt. 

n>l ^° n: \nn<t 

So, making use of the Euler-Maclaurin summation (1.6) to the Cahen's formula we deduce 

r°° He'] 

y2(-l) n - 1 n-"I v (an) = -v / e~ vt T) u {e mu I v (au)) dtdw; 

and repeating the procedure to the second Dirichlet series containing L„(an), the proof is complete. □ 

The next result concerns a hypergeometric integral, which one we integrate by means of Schlomilch 
series of modified Bcsscl and modified Struvc functions. 
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Theorem 8. Let 3ff(V) > and a > 0. Then we have 



/ J v (ax)S(x) — = — f- — r- / — - — -2F 

Jo 1 ' V V 2^Y(u+\)J e at -l 2 



1 1 



t 2 dt 



+ 



ira" [Li 2 (e- a ) + aLii(e- a )] 



2"+ 1 I> + l) 
where Li a (z) stands for the dilogarithm function. 
Proof. Differentiating (f .1) with respect to n, we get 

f 00 2nJ v (ax) dx n(v + l) air 

I ^7 — = r^r- [ivKdn) — L„(on)) — rr \J-,A an ) — L,.(an)) . 

Jo (x 2 +n 2 ) 2 x' y 2n"+ 2 v v ; v 2n v + 1 y uK ' vK " 

Summing up this relation with respect to positive integers n € N, we get 

dx tt{v + 1) 



K v (a) := / J v (ax) S(x) - 



" 2 [I v {an) - L„(an)) 



— n (K( an ) - K(an)) = ^ +2 («) -y^ S# +2 («) • 



n>l 



Now, by the Emersleben-Elbert formula (3.2) we conclude that 



f°° dx 
K v (a) — / Jv(ax) Six) — 
Jo x " 



t 



e* - 1 



J v (ax) sin(xi) 



dx dt. 



Expressing the sine via Ji , we get that the inner-most integral equals 



(3.4) 



f 

Jo 



J u (ax) sin(xi) 



dx = x l- 



nt f°° Jv(ax)Ji(tx) 



dx . 



Now, we shall apply the Weber-Sonin-Schafheitlin formula [63, §13.41] for A = v + \, which reduces to 

if < a < t 



/>OG 

/ J u (ax)Ji(t. 
Jo 



Accordingly, (3.4) becomes 



x)x v 2 dx = < 



2 v +h Vtl> + 1) : 



1 1 



it- . , 



[2^r(, + i) 



if < t < a 



K„(a) = 



/-a j-2 



2-+ 1 I> + i)7 e*-l 



2F1 



1 i 2 



1 1 _ 
2 ' 2 

3 
2 

1 1 



e_\ na v f°° t A 

a 2 ) + 2"+ir(i/ + 1) J a e* - 1 



2^+!r + i) 7 eat - 1 

-,f+2 /•! ^2 



2-f 1 



2 y + 1 r (1/ + i) 7 eat - 1 



■ 2-F1 



t 2 ) dt + 



t + a 



2"+ 1 T(i> + l) J e*+ a -l 



dt 



t 2 dt 



+ ■ 



7ra" 



[Li 2 (e-°)+oLii(e- )] , 



2^+ 1 r(^ + 1) 

where the dilogarithm Li a (z) = J2 n >i znn ~ a , \z\ < 1, has the integral representation 



hi a (z) 



r(a) i 



00 ^a — 1 

e* — z 



dt, 5R(a)>0. 



This completes the proof. 



□ 
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4. Differential equations for Kapteyn and Schlomilch series of modified Bessel and 

modified Struve functions 

Kapteyn series of Bessel functions were introduced by Willcm Kapteyn [32, 33], and were considered 
and discussed in details by Nielsen [41] and Watson [63], who devoted a whole section of his celebrated 
monograph to this theme. Recently, the present authors and Jankov obtained integral representation and 
ordinary differential equations descriptions and related results for real variable Kapteyn series [4, 30]. 

Now, we will consider the Kapteyn series built by modified Bessel functions of the first kind, and 
modified Struve functions 



n>l 



where the parameter space includes positive a > 0, while sequence (a n ) n >i ensures the convergence of 
A^(x). Our first goal is to establish double definite integral representation formula for ^.^(x). In this 
goal we recall the definition of the confluent Fox- Wright generalized hypergeometric function (for 
the general case consult [55, p. 493]): 



(4.1) 



(b,a) 



E 

n=0 



(a) 



pn 



(b)a 



where a, b <E C, p, a > and where, as usual, (A) M denotes the Pochhammer symbol defined, in terms of 
Euler's Gamma function, by 



T(X + p) _ J 1, if /i = 0;AeC\{0} 

T(A) ~ |A(A + l)---(A + n-l), if (j, = neN; AeC 



The defining series in (4.1) converges in the whole complex z-plane when A = a — p + 1 > 0; if A = 0, 
then the series converges for \z\ < V, where V := p~ p a a . 

Theorem 9. Let p > v > and let a e C 2 (R + ) 7 such that a |n = (a n ) ■ Then for 



x e 0, 2min< 1, 



e limsup la^l 1 /" 



— f 



we have 



(2 ' 2) 

Ki) 



xt 



a(s)s^-^ , , 
• 5 S — -, rv dtds . 



Proof. The Sonin-Gubler formula enables us to transform the summands of the Kapteyn series A^(x) 
into 

J vn {xy) 



s^+i E^ 



(y 2 + n 2 ) y" 



dy. 



Making use of the Gegenbauer's integral expression for J a [1, p. 204, Eq. (4.7.5)], after some algebra we 
get 



[ ' vWo VT - ? [ ^! r ("» + 1) Jo v 2 + '■ 



cos(xty) 



dy } dt 



_2_ f 1 1 a n (|(1 -t : 



e -xtn 



n (f(l-t 2 j)_ 



dt 
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where the inner sum is evidently the following Dirichlet series 

* n exp j-n (xt + z/ln ^tzF) ) } 



and p(t) = xt + v\n x ^_ t 2) > for x G (0, 2), since p is increasing on (0, 1). By the Cauchy convergence 
test applied to @ a (t) we deduce that 



C^(l -t^Ye-^ limsupl^l 1 /" < f^LY limsupl^l 1 /" 



< l 



that is, for all x G J a the series converges absolutely and uniformly. By the Cahen formula (1.6) we have 



2> a (t)=\ne a 



d s — -4 r^dzds. 



r (i/a + i) 



Thus 



where the ^-integral 



^ {x) = --= / D s <M*)d*ds, 



has to be evaluated. After indefinite integration, under definite integral, expanding the exponential term 
into Maclaurin series, legitimate termwise integration leads to 



/<M*)<u=(!)"jf(i-< 



2yz-^ e -xzt df 



2r(j/z) 



(!) E 



Consequently 



and thus 



2Y{vz) 
y/n d T (yz + 5) /x 





r f 1 lN i 




) 1*1 


V2' 2) 

{uz,l) 


— xz 



(!) 



•^-J^+iM 



2 T{vz) 

■w 9 i>t + A 



C 1 i) 

V2> 2/ 
(l/Z, 1) 



at v(vt) V2 
(2' 2) 



Ki) 



— ait 



a(s)s ,ys -^- 1 , , 
d s — --r r^dtds. 



r(^ + i) 



The proof is complete. 



□ 



Now, our goal is to establish a second order nonhomogeneous ordinary differential equation which 
particular solution is the above introduced special kind Kapteyn series (1.3). Firstly, we introduce the 
modified Bessel type differential operator 

M[y} = y" + -y'-(l+'^] y; 



this operator is associated with the modified Struve differential equation, reads as follows 

<«» M M s »" + i»'-( 1+ 5)»=;J&)' 
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Theorem 10. Let min(z/, //) > 0. Then for x € J? a the Kapteyn series A — &"(x) is a particular 
solution of the nonhomogeneous linear second order ordinary differential equation 



(4.4) 
where 

*~ (x ) = A r f M ^ rfo + i) /x^ 

""•"^ ; dx J 1 J 1 dt T(vt) \2J 
r fi ^ 

x V2> II 

11 Ki) 



— xt 



a(s)s vs -^- 1 (s-l) , , 

r (i/s + ±) 



Proof. Consider the modified Struve differential equation (4.3) 

M[y] = y"{x) + - y'(x) + y(x) = 



v^I>+i) 

which possesses the solution y(x) — c\I v {x) + c 2 L u (x) + c 3 K u (x) , where K v stands for the modified Bessel 
function of the second kind of order v. Being l v and K v independent particular solutions (the Wronskian 
W\I Vl K v \ = — x^ 1 ) of the homogeneous modified Bessel ordinary differential equation, which appears on 
the left side in (4.3), the choice c 3 = is legitimate. Thus y(x) — I vn (x) — h un (x) is also a particular 
solution of (4.3). Setting v M> vn, we get 

(I vn (x) - L m (x))" + ^ (I vn (x) - li vn (x))' 

/ v \ vn—l 



1 H o~ (tvnix) ~ L vn (x)) 



v 2 j 



Finally, putting x ^ xn multiplying the above display with n~>*a n and summing up in n G N, we obtain 

1 9 a ( xn ) vn 

MiKl = m;w - - («;■(,) - ^ ( ,)) + £ r( J^; m , +1 , 

where all three right-hand side series converge uniformly inside ,f a . Applying the result (4.2) of the 
previous theorem to the series 

KA X ) - &"u+i( x ) = E - (^n{xn) - h un (xn)) , 



n>2 



the summation begins with 2. So, the current lower integration limit in the Euler-Maclaurin summation 
formula related to (4.2) becomes 1. By this we clearify the stated relation (4.4). □ 

In the following we concentrate on the summation of Schlomilch series 

^» : =E^ (Unx)-U(nx)) 

n>l 

■= E ^ nx ) M«*)) • 

n>l 

To unify these procedures, we consider the generalized Schlomilch series like (1.5) 
(4.5) e%(x) = £ ^ (Mxn) - U(xn)) ; 

n>l 

obviously T^'^ (x) , (x) are special cases of &l^(x). However, bearing in mind the asymptotics in 
Sonin-Gubler formula (3.1), we see that the necessary condition for the convergence of &i'^(x) for a 
fixed x > becomes a n — o as n —¥ oo. 
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Theorem 11. Let mm(v , n, x) > and a <E C 1 (M+) be monotone increasing, such that a |n = {a n )n>i, 
and that J2 n >i fl_ '' +l '" la n converges. Then & = &l'^(x) is a particular solution of the nonhomogeneous 
linear second order ordinary differential equation 

X 



m»[6] = m [ 6 ;*] = \ (t*;w)' - 



/+i 



where 



r°° f[e*] 

T a f(x) = nj e - "* y U (a(u)(u' 9 - 1) (^(ra) - L„(xu))) didu. 



Proof. Consider again the modified Struve differential equation (4.3), which possesses the solution y(x) = 
cil 1/ (x) + C2L„(x) + C3K„(x), and choose the particular solution associated with c\ = —C2 = 1 and C3 = 0. 
Transforming (4.3) putting x i-> xn, multiplying it by a n n~^ and summing up the equation with respect 
n € N, we arrive at 



i7l>l / \n>l / 



-y(xn) 



o > rrrVixn) = 

x 2 ^ n^+ 2 1 

n>l 



n>l 



On 



2/ n>\ 



n H—v+l 



Thus 



f 

x 



E 

lTl>2 



~2 



a„(n - 1) 



7 »+i 



y(xn) 



a„(n 2 - f) 



x^ ^ — ' n' i + 2 

n>2 



y(xn) 



(I) 



L/-1 



n u-u+l 



Denote 
(4.6) 



T Q /(x) = £ 



n>2 



y(xn), 



0<(v<^t,x>0. 



Following the same lines of the proof of Theorem 7, by Cahen's formula and the Euler-Maclaurin sum- 
mation we immediately yield the double definite integral representation 

■!>*] 



r°° ri e J 

= H J e - "' y 3„ - 1) y(xu)) dtdu ; 

which immediatelly lead to the stated result. 



□ 



Corollary 1. Let ^ — 1 > v > and x > 0. Then X = Tf/^x) is a particular solution of the nonhomo- 
geneous linear second order ordinary differential equation 



where 



poo Ae 1 ] 

T^(x) = ^y e-"* y 0„ (( U " - 1) - Mxw))) didu. 



Corollary 2. Let \i > v > andx > 0. T/ien 1 = T^(x) is a particular solution of the nonhomogeneous 
linear second order ordinary differential equation 

,2 



M° 



= M 
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where 

T£(sc) = fij e - "' y U (e i7ru (u* 3 - 1)(L„(jcu) - I v (xu))) dtdu. 

Now, a completely different type integral representation formula will be derived for T^ +1 (x) which 
simplifies the nonhomogeneous part of related differential equation (4.7). 

Theorem 12. If u > and x > 0, i/ien we /mue 

Proof. Consider the well-known summation formula [24] 

El 7T 1 

— = — coth(7ra) — —r, a m . 

a 2 + n 2 2a y ' 2a 2 ' r 

n>l 

In conjunction with the Sonin-Gubler formula (f .6) we conclude that 



dt 

7rf 2 / F 



□ 



Remark 1. Actually, the formula 

El 7T f 

— t = — coth(7ra) — a ^ 



a 2 + n 2 2a y ' 2a 2 '' 

n>l 



has been considered by Hamburger [24, p. 130, (C)] in the slightly different form 

00 . _ 00 

(C) 1 + 2V e" 2 ™ a = icot7ria = — + — V — -, a^in; 

y ' ^ Tra it — ' a -\- n 

n—l n—1 

Hamburger proved that the functional equation for the Riemann-zeta function is equivalent to (C), see 
also [11] for connections the above formulae to Eisenstein series. Also, it is worth to mention that further, 
complex analytical generalizations of above formula can be found in [8]. 

5. Novel Bromwich Wagner line integral representation for J u (x) 
As a by-product of Theorem 12, we arrive at the integral relation 

s££+iO»0 = (^( xn ) - M^)) . v > 0,x > 

n>l 12 

which, in the expanded form reads 

r°° ( \\ At f°° 

y M xt) ^ ooth ( 7rt )-_j_ = (|/ + l ) y y e-^ s d u (I u (xu) -U(xu)) dsdu. 
This arises in a Fredholm type convolutional integral equation of the first kind with degenerate kernel 

(5.1) J™ f(xt) (coth(Trf) - 1) = Fv{x) , 
having nonhomogeneous part 

f°° 

(5.2) F v {x) = {v+l) / e-< v+1 > 8 MM*u)-M*u)) dsdu - 

Jo Jo 

Obviously, J„ is a particular solution of this equation. Before we state our result, we say that the functions 
/ and g are orthogonal a. e. with respect to the ordinary Lebesgue measure on the positive half- line when 
Jo°° f( x )9( x ) dx vanishes, writing this as / _L g. 
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Theorem 13. Let v > and x > 0. The first kind Fredholm type convolutional integral equation with 
degenerate kernel (5.1) possesses particular solution f = J v + h, where h <G L 1 (R + ) and 

h(x) _L x-"- 1 ( coth(7ra;) ) , x > 

\ 7TX J 

if and only if the nonhomogeneous part of the integral equation equals F v (x) given by (5.2). 

We mention that h as in the above theorem has been constructed in [16, Example]. To solve the 
integral equation (5.2) we use the Mellin integral transform technique, following some lines of a similar 
procedure used by Drascic-Pogany in [16]. The Mellin transform pair of certain suitable / we define [53] 

/•oo i />c+ioo 

:= / xP-\f(x)dx, jZ-\g);= X->\* p (f)dp, 

Jo Z7T1 Jc-ioo 

where the inverse Mellin transform is given in the form of a line integral with Bromwich-Wagner type 
integration path which begins at c — ioo and terminates at c + ioo. Here the real c belongs to the 
fundamental strip of the inverse Mellin transform 

Theorem 14. Let v > 0,x > 0. Then the following Bromwich-Wagner type line integral representation 
holds true 

/•oo /*[e s 

! dtt 



c+ioo 



/ / e-( v+ V a * u (I v (xu)-L v (xu))dst 
. Jo Jo 



(5 ' 3) Mx) = V " T(*?)r(z* + i)«v- P + 2) - xP ~ Up ' 

where c e (v, u + 1). 

Proof. Applying Jt~ x to the equation (5.1), we get 

Jt v Mxt) ( coth(^) - -L) -^j = J( V {F V ) . 

By the Mellin-convolution property 

■^ P (f*g) = -#v f{rt) ■ g(t) dt^j = Jt p {f) ■ .#i-p(fiO , 

it follows that 

(5.4) Ji v (ar^-^cothTra; - (nx)- 1 )) ■ Jl x _ v {J v ) = Jt v (F v ) . 

The fundamental analytic strip contains (u, 1 + v), because the coth behaves like 

l + i + 0(**), ifz^O 
coth z = < z 3 



l + 2c- 2z [l + C(c- 2z )] if z 



oo 



in both cases 3?(z) > 0. Now, rewriting x v 1 (coth7ra; — {-kx) 1 ) by (C) and using termwise the Beta 
function description, we conclude that 

Jt v {x~ v ~ x (coth the - (the)- 1 )) = - B ( V —^- + l) £( v -p + 2), 

for all p € (v, v + 1). Therefore 



B(E^,^ + 1) C(^-p + 2)' 
which finally results in 

e- {u+1)s U - ~L v {xu)) ds du 



v + 1 f c+[oc ' \Jo Jo 



MX) ~ 2i J c _ ioo B (^,*? + l)<(,/-p + 2) 

where the fundamental strip contains c = v + \. So, the desired integral representation formula is 
established. □ 
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We note that the formula-collection [25] docs not contain (5.3). 



Theorem 15. Let < v < |,a; > 0. Then 



//ere afeo c G v + 1). 

Proof. Consider relation (5.4). Expressing ^i- p {J u ) wo formula [43, p. 93, Eq. 10.1] 

Jt v {J v (ax)) = — - ' , a>0, -v<p<\, 

equality(5.4), by virtue of the Euler's reflection formula becomes 

r(^£ + i)r(^ + i)r(^) c(^-p + 2) 



Jt v {F v ) 



2^r(^ + i) 



2P sin[|(p-^)].r(^ + i)- 

Having in mind that F v {x) is the integral representation of the Schlomilch series % I J v+1 (x), inverting the 
last display by we arrive at the asserted result. □ 
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